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We performed path integral simulations of spin evolution controlled by the Rashba spin-orbit 
interaction in the semiclassical regime for chaotic and regular quantum dots. The spin polarization 
dynamics have been found to be strikingly different from the D'yakonov-Perel' (DP) spin relaxation 
in bulk systems. Also an important distinction have been found between long time spin evolutions 
in classically chaotic and regular systems. In the former case the spin polarization relaxes to zero 
within relaxation time much larger than the DP relaxation, while in the latter case it evolves to a 
time independent residual value. The quantum mechanical analysis of the spin evolution based on 
the exact solution of the Schrodinger equation with Rashba SOI has confirmed the results of the 
classical simulations for the circular dot, which is expected to be valid in general regular systems. In 
contrast, the spin relaxation down to zero in chaotic dots contradicts to what have to be expected 
from quantum mechanics. This signals on importance at long time of the mesoscopic echo effect 
missed in the semiclassical simulations. 

PACS numbers: 72.25.Rb, 72.25.Dc, 73.63.Kv, 03.65.Sq 



I INTRODUCTION 

Spin relaxation in semiconductors is an important 
physical phenomenon being actively studied recently in 
connection with various spintronics applications 0. In 
doped bulk samples and quantum wells (QW) of III-V 
semiconductors at low temperatures spin relaxation is 
mostly due to the DP mechanism 0- This mechanism 
does not involve any inelastic processes, so that the ex- 
ponential decay of the spin polarization is determined 
entirely by the spin-orbit interaction (SOI) and elastic 
scattering of electrons on the impurities. However, in 
case of confined systems such as quantum dots (QD) with 
atomic-like eigenstates, the SOI has been incorporated 
into the structure of the wave functions of the discrete 
energy levels. Without inelastic interactions, an initial 
wave packet with a given spin polarization will evolve in 
time as a coherent superposition of these discrete eigen- 
states. Therefore, the corresponding expectation value 
of the spin polarization will oscillate in time without any 
decay. To obtain a polarization decay in the QD's, ex- 
tra effects have to be introduced into the system, e.g., 
the inelastic interactions between electrons and phonons 
mediated by the spin-orbit 0, 0] and nuclear hyperfme 
effects 0,(1,0. Accordingly, a spin relaxation in QD's 
induced by these effects is a real dephasing process. 

Unlike such an inelastic relaxation in QD's, the DP 
spin relaxation in unbounded systems seems to be a quite 
different phenomenon, because the scattering on impuri- 
ties is elastic and there is no dephasing of the electron 
wave functions in the systems. However, the spin polar- 
ization does decay in time exponentially, as if it would be 
a true dephasing process. To explain this phenomenon, 
let us consider an electron moving diffusively through an 



unbounded system with random clastic scatters. This 
electron is described by a wave packet represented by a 
superposition of continuum eigenstates. During a DP 
relaxation process, the spin expectation value expressed 
as a bilinear combination of these wave amplitudes will 
decay exponentially in time. This process can be easily 
understood from the semiclassical Boltzmann or Fokker- 
Plank approach 2J. Indeed, keeping in mind that the 
SOI has the form cr • h(k), where er is the vector, whose 
components are the three Pauli matrices, and h(k) is the 
effective magnetic field, whose magnitude and direction 
depend on the electron momentum k, one can envision 
spin relaxation as the spin random walk on the surface 
of the unit sphere, similar to that in Fig. 1(c). Starting 
at the north pole, the spin precesses around h(kj) until 
the momentum direction is changed by a scattering on 
an impurity. Thereafter, the magnetic field changes its 
direction to h(k2) and the spin continues its precession 
around this new direction. If the spin rotation angle be- 
tween successive scattering events is small, the sequence 
of such rotations results in a diffusive spreading of the 
initial polarization. 

Returning to QD's, a natural question emerges: what 
sort of spin evolution can be generated by the DP mecha- 
nism in a ballistic QD whose size is much larger than the 
electron wavelength at the Fermi surface and where the 
mean spacing between energy levels is much less than 
fr/T, where T is the mean time between electron colli- 
sions with the boundary? Similar to the example in Fig. 
1, the spin evolution in this semiclassical regime can be 
studied by tracking the spin walk on the sphere, when 
particles move along the classical trajectories inside the 
QD's. Intuitively, one would expect the spin evolution 
in this case to be similar to the spin random walk gov- 
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erned by the impurity scattering in unbounded samples. 
However, this expected analogy with the open system is 
wrong. Indeed, in an unbounded system, the steps of the 
random walk are uncorrelated. This results in a diffu- 
sive decay of the spin polarization down to zero for any 
nonzero SOI. But in case of QD's, the steps of the random 
walk on the sphere are correlated due to the confinement 
of electron trajectories within the dots. As we will show 
below, such correlations not only lead to a spin relax- 
ation much longer than the DP relaxation in unbounded 
systems, but also to a non-zero final polarization value 
at long time for certain quantum dot geometries. Here, 
we do not take into account the inelastic mechanisms 
EJ El IE El which always drive the spin polarization to 
zero in long time. These mechanisms are assumed to be 
absent, because they become inefficient at low sufficiently 
temperatures. 

In this article we carry out a semiclassical analysis of 
the DP relaxation in 2-dimcnsional (2D) QD's of various 
geometries, including a circular dot, a triangular dot, a 
generalized Sinai billiard, and a circular dot with diffu- 
sive scattering on the boundary. We focus on the case 
of the strong SOI, such that the characteristic spin orbit 
length L so = vp /Hh(kF) is not much larger than the dot 
size L. Such a regime can be realized in the InAs based 
heterostructures for L ~ 0.5 — 1/im 7]. We found that 
in the short time scale ~ T the spin relaxation dynamics 
in all geometries shares a common feature: After a fast 
initial drop during the time interval ~ T, the spin polar- 
ization continues to oscillate weakly around some value. 
For weak SOI with L so ^> L, all residual values for differ- 
ent dot geometries are quite close to one up to the cutoff 
time of our numerical simulations (~ 10 3 T). For stronger 
SOI with L so > L, the initial drop of the spin polariza- 
tion is considerably larger compared to the weak SOI 
regime. The spin evolution after that drop depends on 
the dot geometry. In the case of circular and triangular 
dots, which are examples of systems with regular classical 
dynamics, the corresponding spin polarizations approach 
nonzero residual values. However, in the case of chaotic 
and random systems (e.g., Sinai billiard and circular dot 
with rough boundaries, respectively), the spin polariza- 
tions slowly decrease to zero after that initial drop. But 
this decreasing is much longer than the DP relaxation in 
an unbounded system, in which the mean impurity scat- 
tering time is ~ T. For very strong SOI with L so < L, 
the spin polarization after the initial drop reaches zero 
and later on oscillates with a large amplitude. 

These results clearly demonstrate that the spin evolu- 
tion in QD's is qualitatively distinct from the DP spin re- 
laxation in unbounded systems. In order to elucidate the 
physical origin of this phenomenon, two investigations 
have been performed. First, the spin evolution along a 
single electron trajectory was studied in detail, which 
provided a clue for understanding the above-mentioned 
polarization behavior. Second, the residual polariza- 



tion obtained from the classical simulations for a circular 
quantum dot was compared with that derived from the 
exact solution of the Schrodinger equation. A good agree- 
ment between the results from these two approaches has 
been found. However, for QD's with chaotic and ran- 
dom electron dynamics, the general quantum mechanical 
analysis revealed a contradiction to the long time spin 
evolution observed in our semiclassical simulations. 

The article is organized in the following way. In Sec- 
tion II the general expression of the polarization will be 
derived for the spin evolution via classical path integrals. 
In Section III the results of the numerical simulations in 
different quantum dots will be demonstrated. The quan- 
tum mechanical theory for the spin polarization in the 
circular quantum dot will be presented in Section IV, 
with the calculation in detail shown in the Appendix. 
Discussion and conclusion will be given in Section V. 

II PATH INTEGRALS FOR THE SPIN 
EVOLUTION 

The Hamiltonian of the system, 

H = F + T-h(k), (1) 

consists of the spin independent part Ho, which is the 
electron kinetic energy plus the 2D confining potential 
V(r), and the spin-orbit interaction. In III-V semi- 
conductor heterostructures the effective "magnetic" field 
h(k) is given by the sum of the Rashba || and the Dres- 
selhaus |9| terms. If the z-axis is chosen perpendicu- 
lar to the heterostructure interface, the magnetic field 
Iir contributing to the Rashba term has two compo- 
nents (h R (k) , h v R (k)) = (aRky, — ctjik x ), where %k = 
(hk x ,hk y ) is the momentum operator. In the 2D con- 
finement, the magnetic field ho contributing to the Dres- 
selhaus term contains both linear and cubic parts with 
respect to k 0. In a [001] oriented QW the linear term 
has the components {h x D {k), H u D (k)) = (aokx^—anky). 
For heterostructures with a typical ~10nm confinement 
in z-direction, the linear part of ho is usually larger than 
the cubic part, except the case of high doping concen- 
tration [ll|. The Rashba term is not zero only in het- 
erostructures with asymmetry in their growth direction. 
This term can be much larger than the Dresselhaus term 
in the narrow gap InAs based systems ■ In this article 
we will study the spin evolution induced by the Rashba 
term. But, since the SOI Hamiltonians corresponding 
to the Rashba and the linear Dresselhaus terms can be 
transformed from one to the other by the unitary matrix 
(<j x + a v )/y/2, our results are also valid for systems in 
which the linear Dresselhaus term dominates the SOI. 

Let us suppose E n to be the n-th quantized energy level 
with the eigenfunction ip n , which is a two component 
spinor. At zero magnetic field this quantum state is at 
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least doubly degenerate. Let 

ip(r) = e ikr $(r - R) X 



(2) 



be the wave packet created at time t = 0, centered at 
the point R, and propagating with the 2D wavevector 
k. The function <!>(r) is assumed to be slowly varying 
within the scale of the electron wavelength 27r/fc and 
normalized, so that the integral J |$(r)| 2 d 2 r over the 
QD volume is equal to 1. The initial spin polarization 
P(0) = Xa^afiXp is the sum over the two compo- 
nents Xa of the spinor %, where a e {1,2}. For t > 0, 
the wave packet evolves in time as 



V>(r, t) 



c n ip n (r)e- lE - t/h 



where 



V?n(rMr) d 2 r. 



(3) 



(4) 



In terms of "0(r, t) the time dependent spin polarization 
is expressed as 

P = £ / ^( r ' QvbpM** *) d2r > ( 5 ) 



a/3 



with three components P(i) = {P x {t),py{t),P z (t)). 

For further analysis it is convenient to introduce the 
retarded and advanced Green functions, 



G r a Jt-t',r,r')=GZ(t'-t,r',r 



p K u - u ,1,1 j — u^^i - fc, i ,i; (6) 



which are 2x2 matrices acting on the SU(2) spin space, 
where 0(i — t') is the Heaviside function. Using these 
Green functions, the spin-spin correlation function can 
be defined as 

7f«(r,r';(-i') (7) 

where i,j € {x,y, z}. This definition together with Eqs. 
(J3EJ| lead to the expression for the polarization evolution 
in time, 

P i (t) = - j i^ (r,r';i)$(r-R)<r(r'-R) 

xe ik ( r - r ')pi(o)d 2 rdV. (8) 

For classical simulations below, the semiclassical ap- 
proximation of Eq. © is required. It can be derived 
from a standard path integral formalism [l2T |. by repre- 
senting the retarded Green function in Eq. Q as the 



sum of products, 
G r ag (t-t',r,r') 



(9) 



dri • • • dr„ 



E 



(r 1 ,ai|e-^ ( * 1 -* 2 )|r2,a 2 ) 



-iH(t n -t'),i 



of the evolution operators e -1 -^** - *" within the infinites- 
imally short time intervals (tj — tj). Thereafter, the 
Green function can be expressed as the path integral of 
T exp [j-S(t — t',r, r')l , where the action 



S(t-t',r,r') 



-v 2 (r)-V(v(r)) 



m*v(r) 



(10) 



is a time integral of the particle Lagrangian evaluated 
along a trajectory starting from r' at time t' and ending 
with r at time t, where v(t) = In this Lagrangian, 
the constant term m*a 2 R j2 is ignored, because it only 
gives a phase factor. Since the SOI Lagrangians on dif- 
ferent parts of the trajectory do not commute, one has 
to keep different exp \i S(t — t', r, r')] in the order of the 
sequence in Eq. which is preserved by the time or- 
dering operator T. 

By using the saddle point approximation, the path in- 
tegral in Eq. (Q can be reduced to a sum over all classical 
trajectories 7 jl2|, 

G r (t-t',r,r') (11) 
= ^J2 ^fJM^ So{t - t '^ T Mt - t', r, r') , 

Z7T * — ' 
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with the spin independent monodromy matrix J(r, r') = 

det ( gl f 0°i ^ an d the spin independent classical action 

So(t — t',r,r') along the classical trajectories. The spin 
dependence part of the Green function is represented by 
the unitary matrix 



U(t 



= Te-*/>(^)- 



'dr. 



(12) 



Such a decoupling of the spatial and spin degrees of free- 
dom can be done under the assumption that the classical 
paths are only weakly perturbed by SOI, which is rea- 
sonable, when the SOI parameter cxr is much less than 
the electron Fermi velocity. Under this assumption, all 
quantities J, S, and U are evaluated on the unperturbed 
trajectories. 

Inserting Eq. ifTTjl and fL^I into Eq. J2J) and JBJ), we ob- 
tain a semiclassical expression for the spin polarization. 
This expression can be substantially simplified after in- 
tegrating over coordinates r and r' in Eq. JSJ. Indeed, 
let us consider the integral in Eq. (JHJ, 

J y/J(r", r)e^ So( *' r "' r) $(r - R)e lkr [/(i, r", v)d 2 r. (13) 
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In the semiclassical limit, the exponential function 
exp [jrSo(t, r",r)l rapidly oscillates as a function of r 
with a period given by the Fermi wavelength. However, 
J, U, and $ are slowly varying functions of r. The length 
scale of J's variation is given by the dot size. The spa- 
tial changes of U are controlled by the spin orbit length 
L so = h/(m*aii), which is assumed to be much larger 
than the Fermi wavelength. Therefore, the influence of 
the SOI on the saddle-point position can be ignored. The 
variation of $ also can be ignored, because this function 
was assumed to change weakly within the length scale 
equal to the electron wavelength. Under these approxi- 
mations, we obtain the saddle-point equation in the form 



dS (t,r",r) 
dr 



hk = 0. 



(14) 



This equation is the classical equation of motion. It de- 
termines the trajectory r = r (r"(i), p(0)) which passes 
through the given point r"(t) at the instant t, on condi- 
tion that at t = the initial momentum was p(0) = ?ik. 
Therefore, the saddle point r is a particle coordinate at 
t = belonging to this trajectory. Since the integral over 
r' in Eq. JHJ is taken around this extremum, the value 
r' = r = ro are inserted into all slowly varying functions 
J, U and 

Further, to calculate the integral over r in Eq. (|13fl . 
the action 5o(£, r",r) is expanded around r = r up to 
the second order, 

S (t,r",r) + hk = S (t, r",r ) 
£ or l or J 

The integration over r and r' in Eq. (JSJ gives 
(27r) 2 /det ( 85 ° (t 'f f o) Y Combining this Jacobian with 
J(r", r ) we obtain 



det 



aSo(*,r",r ) 
dr" l dr 3 



det 



dS (t,r",T ) 
dr l dr 3 



det ' 9r ° 



dr"3 



By using the identity 



(16) 



(17) 



Eq. Q can be integrated over r , instead of r", which 
leads to the expression of the semiclassical spin polariza- 
tion, 

P'cit) = J ^'(r,r',t)|$(r'-R)|W, (18) 



Equation i|18|l describes the spin evolution of a particle 
initially distributed around the point R with the proba- 
bility density |$(r'— R)| 2 . This particle starts its classical 
motion from the point r' with the momentum Kk at time 
zero and arrives in the position r at time t. In the fol- 
lowing, we are interested in the spin evolution averaged 
over an ensemble of electrons with uniformly distributed 
coordinates R and random directions of the initial mo- 
menta on the Fermi surface. After averaging Eq. Ijl8(l 
over R and the angular coordinate (9k of the momentum 
k, we obtain the simple expression: 



Kit) 



P J (0) 
An 



J R ij (r,r',t)d 2 r'd6 k . 



(20) 



It should be noted that after the integration over R this 
expression does not depend on the initial wave packet 
envelope $(r — R). Therefore, the same Eq. Ij20|l holds 
for $ = const, so that the initial state can be simply a 
plane wave. 



Ill NUMERICAL RESULTS 

Equation (|2UI) is the basic equation for our numerical 
simulations of the spin polarization. Below we will re- 
strict ourselves to the case when the initial polarization 
P(0) is directed along the z axis, so that P 2 (0) = 1, 
and the polarization to be calculated at time t is also in 
^-direction. 



III.l Spin evolution in ballistic quantum dots 

Consider a free electron confined inside a quantum dot 
and moving along the trajectory 7, which consists of the 
successive straight segments 7j of the lengths lj with j = 
1, 2, 71. The spin state along this trajectory can be 
described by the evolution operator U 1 — U (t, r, r') in 
Eq. (|12|l with t' = 0. This operator can be represented 
as a product 



of the individual operators 



U- 



exp [-itpj Jj] . 



(21) 



(22) 



with ipj — 



Jj = Nj • a. Thereby, N, = rij x e z 



with 



is the unit vector parallel to the effective magnetic field 
h(k) = Q!ij(k x e z ), where rij — k/|k| is the unit vector 
along the trajectory segment j and e z is the unit vector 
in ^-direction. Since Jj is a vector in the space of the 
Pauli matrices, the individual operator in Eq. I]22ll has a 
simple form 



R l3 (r,r',t) = Tr [ CT ^(t, r, rV't/^, r, r')] . (19) 



U- 



cos(i/>j) 1 — zsin(^j) Jj, 



(23) 
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with the identity matrix 1. 

Let us assume the j-th segment jj to have the angle 
Wj with respect to the x-axis. Accordingly, the vector 
Nj has the angle Wj — tt/2, so that we get the explicit 
expression Jj — sxn(wj)a x — cos(wj)a y . In SU(2) rep- 
resentation, the operator U lj can be expressed as the 
matrix 



cos(V'j 



-sin(^)e l w i 
which acts on the spin state 



sin(V>j)e 1 
cos(ipj) 



X 



Xi 

X2 



cos(6l/2)e^ 1 
sin(6»/2)e i ^ 



(24) 



(25) 



In SO (3) representation, the operator U lj corresponds 
to a spin rotation around the axis Nj through the angle 
2ipj . The three components of the spin expectation value 
are related to the spinor x by 



2Re(x lX2 ) 
2Im(x!x 2 ) 



(26) 



1X2 1 



For convenience, we will call the vector projections s, S 
[—1, 1] as spin components, although they are twice 
larger than the corresponding values for the spin 1/2. 

As an example of spin evolution induced by the Rashba 
interaction, let us consider an electron confined inside a 
quantum dot in Fig. 1(a), moving along the trajectory 
7 which consists of three straight segments 71, 72, and 
73 with the respective lengths l\, I2, I3 and the angles 
u>i = tt/2, u>2 — 7T, u>3 = 37r/2. The initial spin state 
of this electron is polarized in ^-direction, which is rep- 
resented by an arrow in Fig. 1(c). This arrow is pro- 
jected down to the origin (0, 0) on the s x s y plane in Fig. 
1(b). When the electron starts its motion from the initial 
point p along the segment 71 (Fig. 1(a)), its spin rotates 
around the axis Ni = (1,0,0) and circumscribes an arc 
on the 3-dimensional sphere in Fig. 1(c). This curve is 
projected down onto a straight line on the s x s y plane. 
This line is parallel to 71, but runs in a direction oppo- 
site to 71, as shown in Fig. 1(b). After the first collision 
with the boundary the electron further moves along the 
segment 72, while its spin rotates around N 2 = (0,1,0) 
and circumscribes the second arc on the sphere in Fig. 
1(c). The spin projection in Fig. 1(b) now runs parallel 
to 72 in the direction opposite to electron motion along 
72. It is easy to see that the spin evolution on other seg- 
ments follows the same rule: When an electron passes 
through the j-th segment in a certain direction, the spin 
circumscribes on the 3D unit sphere an arc around the 
axis Nj. This arc, in its turn, is projected onto the s x s v 
plane as a straight line parallel to the electron trajectory, 
but oppositely directed to it. 

Further, let us proceed from the spin evolution on in- 
dividual trajectories to the spin evolution averaged over 
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FIG. 1: (Color online) (a) Electron motion inside a quan- 
tum dot. The trajectory consists of three straight segments 
71, 72, and 73. (b) The corresponding spin evolution on the 
s x Sy plane, which is projected from (c). (c) The spin evo- 
lution induced by the Rashba spin-orbit interaction on the 
3-dimensional unit sphere. 



an ensemble of trajectories. We consider an ensemble of 
electrons distributed uniformly within a bounded area of 
a 2-dimensional heterostructure. At t = these electrons 
have random outgoing angles but the same spins polar- 
ized in z-direction. Let si (t) be the z component of the 
electron spin at time t for the i-th trajectory. Then, in 
our numerical simulations the integral in Eq. i|20|) can be 
replaced by the sum, 



1 ™ 



(*), 



(27) 



where the sum runs over n individual trajectories. The 
so averaged spin polarization will be calculated in the 
following five systems: 

(a) In 2-dimcnsional bulk (Fig. 2(a)) with the elastic 
collision length I distributed according to the Pois- 
son law Prob(Z) = e~ l / lm /l m , where l m is the mean 
free path. It is a stochastic open system. This is 
just the system where the conventional D'yakonov- 
Perel' spin relaxation has to be observed. 
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(b) In a ballistic circular quantum dot of radius 1 with 

the smooth boundary in Fig. 2(b). Since the 
boundary is smooth, the incident and reflection an- 
gles on the boundary are the same. Since the sys- 
tem is ballistic, no scattering occurs inside the dot. 
It is an integrable system with a high spatial sym- 
metry. 

(c) In a ballistic triangular quantum dot with the 

smooth boundary in Fig. 2(c). It is an integrable 
system of lower symmetry compared to the circular 
dot. 

(d) In a generalized Sinai billiard with the smooth 

boundary in Fig. 2(d). It is a deterministic but 
strongly chaotic system. The boundary geometry 
generates an ergodic dynamics in the phase space. 

(e) In a ballistic circular quantum dot like Fig. 2(b), but 

with random reflections from the boundary. The 
reflection angle takes random values between ~n/2 
and 7r/2 with respect to the boundary normal. It 
is a stochastic closed system and corresponds to a 
quantum dot whose boundary is not perfect in the 
scale of the electron Fermi wavelength. 

The mean free path l m in bulk in Fig. 2(a) is set to 1. 
The sizes of the triangular and Sinai dots, as shown in 
Fig. 2(c) and (d), are chosen to be \/2tt ss 2.5066 and 
y/327t/(16 — 7r) « 2.7961, such that these dots have the 
same area ir as that of the circular dot in Fig. 2(b). We 
will use the dimensionless time unit, such that during the 
time interval 1 a particle moving with the Fermi velocity 
travels a distance of the length 1. 

III. 2 Results of the numerical simulations 

In Fig. 3 the time dependences of P* (t) for 2124 elec- 
trons in the open system (Fig. 2(a)) with L so = 10, 6 and 
2 are plotted by solid curves Ci, C*2, and C3. One can 
see that the relaxation time increases with L so . These 
curves can be fitted by the well known expression for the 
longitudinal DP relaxation 13], 

P D „(t)=exp(=i^V (28) 

which is shown by the dashed curves in Fig. 3. This ex- 
pression was derived under the assumption of sufficiently 
large L so 3> l m . For not so large L so the fitting is not 
good, as it can be seen for the curve C3 around its first 
drop at t = 4. In this regime the spin rotates rather 
fast, so that most of the spins si\t) evolve to negative 
values before the electrons encounter their first collisions 
with impurities. Therefore, P c z (i) can evolve to a deep 
negative value within a short time interval. But later on 
PI (t) approaches to the asymptotic value P c z = (curve 




-2 2 _10 1 




1 2.5 1.4 2.8 



FIG. 2: (Color online) Electrons trajectories (solid lines) for 
short time intervals: (a) in bulk, (b) circular quantum dot, 
(c) triangular quantum dot, and (d) Sinai quantum dot. 

C3 in Fig. 3). These results from Monte Carlo simula- 
tions confirm the well known DP relaxation in unbounded 
systems. 




' 2 4 20 40 60 80 100 



Time t 

FIG. 3: Solid curves Ci, C2, and C3 represent the time depen- 
dent polarization Pc(t) for 2124 particles in an unbounded 
QW with L ao = 10, 6, 2 and the mean free path l m = 1. 
The particles were initially placed inside a circular area of 
the radius R — 1 and polarized in ^-direction. The dashed 
curves depict the DP relaxation calculated from Eq. 11281 . For 
comparison, curve C4 shows P<f (t) for 2124 particles confined 
inside a circular dot of the radius R = 1 and L ao = 10. 

If electrons are confined inside the smooth circular dot 
(Fig. 2(b)), the relaxation of P*{t) is considerably sup- 
pressed, so that at large L so the spin polarization remains 
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close to 1 at large times, as the curve C4 in Fig. 3 demon- 
strates for the case of L so = 10. At this regime, the 
suppression of relaxation takes place in all other quan- 
tum dots, like the circular dot with the rough boundary 
(curve Cq), the triangular dot (curve C7), and the Sinai 
billiard (curve C$) in Fig. 4. In all of these curves the 
P^{t) values fall into the range between 0.97 and 0.98 at 
large times up to t = 10 3 . 

On the other hand, the spin polarization evolves very 
fast down to if L so is smaller than the dot size. The 
corresponding time dependence of P* (i) is similar to that 
shown in Fig. 3 (curve C3), with a sharp drop at the 
beginning followed by oscillations around zero. 




20 40 60 80 100 



Time t 

FIG. 4: (Color online) Time dependence of P*(t) for L ao = 10 
in the smooth circular dot (curve C5), the circular dot with 
the rough boundary (curve C§), the triangular dot (curve CV), 
and the Sinai billiard (curve C%). 




Time t 

FIG. 5: (Color online) Time dependence of Pf (t) for L so — 2 
in the smooth circular dot (curve C5), the circular dot with 
the rough boundary (curve Cg), the triangular dot (curve CV), 
and the Sinai billiard (curve Cg). 



For an intermediate L so the spin relaxes according to 
different scenarios, depending on the quantum dot geom- 
etry. As an example, Fig. 5 shows the function P*(t) for 
various dot geometries at L so = 2. After a fast initial 
drop, the polarization further relaxes to in the Sinai 
billiard (curve Cg) and in the circular dot with the rough 
boundary (curve C§). However, in the smooth circu- 
lar (curve C5) and triangular (curve CV) dots this func- 
tion oscillates around a constant value at large times. It 
should be noted that in the former two examples the spin 
polarization relaxes to zero at much longer times than the 
DP relaxation time in the unbounded system (Fig. 3), 
although the mean elastic scattering length there is com- 
parable to the dot size. The relaxation times for C§ and 
Cg, in Fig. 5 increase rapidly with higher L so . Thus, at 
L so = 10 we could not detect any systematic decrease of 
the spin polarization in the Sinai billiard and rough circu- 
lar dot, up to t — 10 3 , which is by an order of magnitude 
larger than the range plotted in Fig. 4. 

An interesting feature of P c z (i) in the regular systems, 
like the triangular and smooth circular dots, is the appar- 
ent oscillation of the polarization. It can be seen at Figs. 
4 and 5, although the oscillations in the latter figure are 
more profound for the case of the circular dot, compared 
to almost vanishing ripples in the triangle. These oscilla- 
tions do not disappear at large times and their amplitude 
increase with the strength of SOI. We can not say much 
about their nature. Probably, they are associated with 
the role of periodic trajectories in regular systems. A 
special study is required to understand the origin and 
characteristics of these oscillations. 



Q. 




L 

so 

FIG. 6: The residual polarization P c z vs the spin rotation 
length Lso for a smooth circular dot. 

At long time the spin polarizations in both regular 
quantum dots (triangle and smooth circle) in Fig. 4 and 
5 oscillate around certain nonzero residual values. These 
residual polarizations P c 2 are L so dependent, as plotted 
in Fig. 6 for the circular dot. 
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III. 3 Spin evolution along individual trajectories 

The existence of the nonzero residual polarization in 
regular quantum dots and long spin relaxation time in 
chaotic systems are fundamentally distinct from the DP 
spin relaxation in the boundless QW. Such a distinction 
is surprising, because at first sight the spin walks on the 
sphere in Fig. 1(c) should be randomized by scattering of 
particles from dot boundaries, similar to randomization 
by impurity scattering in unbounded systems. However, 
this simple point of view is wrong, because there is an 
important difference between the impurity scattering and 
the boundary scattering. For convenience, let us define 
the scattering with a direction change smaller than ir/2 
as a 'forward' scattering and that larger than ir/2 as a 
'backward' scattering. If the particles are isotropically 
scattered by an impurity, half of them continue to move 
'forward'. However, if the particles are scattered by a 
smooth boundary, the particles with incident angles be- 
tween — 7r/4 to 7r/4 with respect to the boundary normal 
will be reflected 'backward'. Since statistically more par- 
ticles hit the boundary within this range of angles, the 
'backward' scattering prevails in DQ's. This property 
of particle scattering can also be extended to QD's with 
rough boundaries. Further, according to Fig. 1, a 'back- 
ward' particle motion is mapped onto a 'backward' spin 
walk. Hence, if the spin moves away from the north pole 
in Fig. 1, after a boundary scattering the spin is more 
likely bounced back towards the north pole. Such a non- 
Markovian statistics of the spin walks gives a clue for 
understanding the numerical results in subsection III. 2. 

In order to make this argument more clear it is instruc- 
tive to study in detail the spin evolution along a single 
trajectory. As described in Fig. 1, the spin motion on the 
unit sphere can be projected onto the s x s y plane. After a 
long time the spin path on the sphere will cover a region 
and produce a certain pattern on the s x s y plane. In the 
circular dot this pattern looks rather ordered. If the elec- 
tron moves along a triangular periodic trajectory (Fig. 
7(a)), the pattern is a rounded triangle (Fig. 7(i)). If the 
trajectories are hexahedral and star-like (Fig. 7(b) and 
(c)), the corresponding patterns are a rounded hexagon 
and a rounded star (Fig. 7(j) and (k)). If the trajecto- 
ries are non-periodic, e.g., Fig. 7(d), the pattern is a disc 
(Fig. 7(1)). A common feature of these patterns is that 
they have the same size, which is less than 1 in the case 
of L so — 5. These patterns are highly stable up to the 
observation time t = 10 4 . It implies that the spin on the 
unit sphere cannot move far away from the north pole, 
so that s^\t) cannot take negative values. Our analysis 
of various trajectories with various initial conditions has 
confirmed this general feature of the spin evolution in the 
circular dot. Hence, a non-zero P c z in Fig. 6 at infinitely 
long time is obviously expected. 

In the triangular dot, two periodic and one non- 



periodic trajectories are shown in Fig. 7(f), (g), and (h). 
The corresponding spin patterns (Fig. 7(n), (o), and 
(p)) are less symmetric and have less predictable sizes 
than those in the circular dot. For the trajectory in Fig. 
7(g) the pattern in Fig. 7(o) touches the circular border. 
Nevertheless, our investigation shows that the patterns 
of most other trajectories are quite stable up to the ob- 
servation time 10 4 and do not touch the boarder. On 
this reason the spin polarization being averaged over tra- 
jectories is expected to relax to a positive residual value, 
although this value is smaller than that in the smooth 
circular dot. 




-1x1 -1 S 1 x 2.5 -1 S 1 



FIG. 7: Electron trajectories on the xy plane ((a)-(h)) and 
respective spin evolution patterns on the s x s y plane ((i)-(p)) 
for Lso = 5. (a), (b), and (c) Periodic triangular, hexahedral, 
and star-like trajectories in the smooth circular dot. (d) A 
non-periodic trajectory in the smooth circular dot. (e) A 
stochastic trajectory in the rough circular dot. (f) and (g) 
Two periodic trajectories in the triangular dot. (h) A non- 
periodic trajectory in the triangular dot. 

In the circular dot with a rough boundary, the reflec- 
tion angles are stochastic, as shown in Fig. 7(e). Within 
the observation time t = 10 3 the corresponding spin pat- 
tern on the s x s v plane has spread out to a much larger 
area (Fig. 7(m)) than those in the smooth circular dot 
(Fig. 7(i)-(l)). Furthermore, the pattern in (Fig. 7(m)) 
is still expanding. The corresponding spin state on the 
3D sphere can penetrate into the lower hemisphere after 
t = 10 3 . However, it can return back to the north sphere 
again. Therefore, the z component of this spin state os- 
cillates between negative and positive values. When av- 
eraged over many trajectories, such oscillations sum up 
to a relaxation curve, similar to Cq in Fig. 5. 

In the Sinai billiard, the s x s y pattern resembles that 
in the rough circular dot. Consequently, the spin relax- 
ation dynamics in both cases have similar characteristics 
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(curves C& and C$ in Fig. 5). 

A general trend seen from Fig. 7 is that the confine- 
ment of the particle motion in QD's makes the spin to 
be also confined within the upper hemisphere, if L so is 
larger than the size of the QD's. For a smooth circular 
dot, this trend can be easily understood from the 'back- 
ward' scattering effect described at the beginning of this 
subsection. Since all trajectories in this case have a sim- 
ple geometry, one can easily see that particles are more 
frequently scattered from the boundary in a 'backward' 
direction. But although this argument holds for general 
bounded systems, it is less evident for other QD's besides 
the smooth circular dot. In general case, the trend to- 
ward the spin confinement can be argued in a different 
way: As seen from Fig. 1, the projected spin path on 
the s x s y plane in Fig. 1(b) is more or less a rescaled 
curve of its particle trajectory in Fig. 1(a). But in real- 
ity the mapping from a trajectory to the corresponding 
spin path is not simply a rescaling, because the spin rota- 
tions on the sphere are non-commutative. For example, 
a closed particle trajectory is in general mapped onto an 
open spin path. However, if L so is large, the spin path 
is restricted to a small part of the sphere. According 
to Eq. H21I22JI . a closed particle trajectory produces an 
open spin path of the linear size ~ l/L so , while the dis- 
tance between the initial and the end points of the path is 
only ~ l/ig - The mapping between the trajectories and 
the spin paths is then similar to a mapping between two 
Euclidean spaces. Therefore, with the accuracy l/L 2 so , 
the spin paths are the rescaled particle trajectories and 
those paths are confined because the particle trajectories 
are confined. It should be noted that such a tendency 
for the spin confinement turns out to be strong even for 
not so large L SOl as one can see from the spin dynamics 
shown in Fig. 5 for L so = 2. 

The above argument about the spin confinement does 
not take into account a long time evolution. Even at 
large L so , small corrections due to non-commutativity 
of spin walks will accumulate in time. As a result, the 
spin can slowly drift toward the lower hemisphere. The 
expanding pattern in Fig. 7(m) of the rough circular dot 
is an example of such a long time behavior. However, 
in contrast to that unstable pattern, the patterns from 
regular systems (Fig. 7(i)-(p) besides (m)) remain stable 
in time. This difference between the single trajectories of 
random and regular systems is consistent with the spin 
relaxation curves shown in Fig. 5. 

Such a distinction between regular and chaotic sys- 
tems follows from fundamental properties of regular and 
chaotic systems. It can be understood from consideration 
of periodic orbits. After a particle runs along a periodic 
orbit 7 and completes a period, its initial spin state x 
will evolve to U 7 x with U~ i = exp[— iilRer], which repre- 
sents a rotation around the axis R through the angle 20. 
Both R and Q are determined entirely by the geometry 
of 7 and by the value of L ao . After the particle repeats w 



periods, all spin positions (U.y) w x, corresponding to the 
end points of all periods w — 1, 2, • • -, are located on a 
closed circle. This circle can be obtained by rotating the 
north pole around R, if the initial x is related to the spin 
polarized in the north pole direction. The other points 
on the periodic orbit are mapped onto spin states around 
this circle. Taking many periodic orbits into account, one 
obtains a set of different axes R and consequently a set of 
circles passing through the north pole. Hence, when av- 
eraged over all periodic orbits, spin spends more time in 
the upper hemisphere. This means that at least the fam- 
ily of the periodic orbits contributes to a nonzero residual 
polarization. How significant is this contribution to the 
whole residual value depends on the amount of the peri- 
odic orbits in a system, which is quite different in regular 
and chaotic systems. In a regular system the family of 
periodic orbits has a finitely positive measure and a bun- 
dle of adjacent nearly periodic orbits. These adjacent 
trajectories behave like periodic orbits if the time is not 
too large, because their linear deviation in time from the 
periodic orbits is small. On the contrary, the periodic 
orbits in chaotic systems are of zero measure 14]. Fur- 
thermore, their adjacent trajectories deviate from them 
exponentially fast. Therefore, with increasing time, the 
weight of the periodic orbits and their adjacent trajec- 
tories becomes exponentially small in chaotic systems, 
while it is a nonzero value in regular systems. Hence, 
as long as we consider only periodic orbits, the residual 
spin polarization has to be a positive number for regular 
systems and zero for chaotic systems. 

The individual trajectory study in a larger time scale 
carried out in this subsection helps us to understand 
some of the results in subsection III. 2. However, al- 
though the existence of the nonzero residual polarization 
P c z is apparent from Figs. 3-6, one cannot exclude a 
possibility that P c z will decay to zero in a much larger 
time scale, since the numerical simulations in all these 
figures are truncated within a finite time. Therefore, we 
can not definitely answer the question whether the spin 
polarization relaxes to zero at the infinitely long time, 
or to a nonzero residual value. For the smooth circular 
dot the latter alternative is corroborated by an analysis 
of the spin polarization from the exact solution of the 
Schrodinger equation, as shown in the next section. 



IV QUANTUM MECHANICAL POLARIZATION 
IN THE CIRCULAR QUANTUM DOT 

Due to the time reversal symmetry, the quantized en- 
ergy levels E n of the Hamiltonian H in Eq. 0) are, at 
least, two fold degenerate with the corresponding spinor 
eigenfunctions <£> n a; where a €E {±} is the degeneracy in- 
dex. In the basis of these states a normalized wave func- 
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tion il){v,t) can be expanded as 



^r,i) = ^c na ^ la (r)e- 4B "*/\ (29) 



with the coefficient 

c«a = / vLWV'W d 2 r . 



(30) 



The expression of ^(r, t) in Eq. <|29(l differs from Eq. © 
only by the degeneracy index a, which is explicitly writ- 
ten here for convenience of our further analysis. Taking 
the notation 



Vw(r, t) = c na ^„ a (r) e i Ent/h , 



(31) 



and V'na(r) = -0 na (r, 0), the z component of the quantum 
mechanical polarization in Eq. J^J can be expressed as 

P z (t) = {^(v,t)\a z \^{v,t)) 

= E/^a(r)a>„ b (r)d 2 r (32) 

nab 

71^771, ab 

The first sum in this equation is time independent, 
while the second sum oscillates in time, so that its av- 
erage over a sufficiently long time interval turns to zero. 
It is interesting to find out whether the former term co- 
incides with the residual polarization in Fig. 6. Such 
a coincidence is not evident because the time dependent 
sum can give rise to large variations of P z {t) after long 
time t. Moreover, the semiclassical theory employed in 
the previous section can be not valid for times larger than 
the mean distance between energy levels near the Fermi 
energy. We can check such a coincidence at least for the 
simple case of a circular dot with the smooth boundary, 
by calculating the residual polarization 



(33) 



because the analytic solution of the Schrodinger equation 
with the arbitrarily strong Rashba interaction is available 
flif . In this section only the key steps of the calculation 
are presented, while the calculation in detail is shown in 
the Appendix. 

Let us consider a circular quantum dot of radius R 
with the confining potential 



V{p) = 



for < p < R 
oo for R > p 



(34) 



written as a function of the polar coordinates r = r 
The eigenfunctions of the n-th eigenvalue E n are [l 



¥>n+( r ) 



Up) 
s ^ +1 >W(p) 



(35) 



and 

<Pn-( T ) 
where the function 

9 AO 



-e-*"*/*(p) 



-au JAKO + JAk-0 

a v -iJu(k+0 + Ju(k-0 



(36) 



(37) 



contains the z^-th order Bessel functions of the first kind 
</„(£), the normalization constant d, the parameters 



J v {k~) _ J v+ x{kJ) 
J»(k+) J v+1 (k + )' 



the wave numbers 



+AsTb 



(38) 



(39) 



and the index 

v = A - 3/2 with A = 1, 2, • • • . (40) 

Therein, the dimensionless parameters £ = p/R, £ — 
2m*ER 2 /h 2 , and b — 2a.Rm* R/Ti 2 have been used. The 
wave numbers k± are quantized because the energy levels 
e are determined by the zeros of the function 

Z v {£) := J u (k-)J v+1 (k+) + J v (k + )J v+1 (k-). (41) 

We chose the plane wave 

lf>(r) = ( I ) e lkr (42) 

as the initial state. After inserting <p n+ (r) from Eq. I|35|l 
and (p n -(r) from Eq. (|36|) together with Eq. I|31|) into 
Eq. i|33|) and averaging over directions of the vector k we 
obtain 



P z = 2nJ2 (|cn+| 2 - \c n A 2 ) {Fu - G n ), 



(43) 



with 



F„ = d 2 

G n = d 2 



2 r(l) 



ail 



r(2) 



r(3) 



(44) 



where the coefficients ii , li 2 \ and i£ are presented 
in Eq. The coefficients |c„±| 2 in Eq. can be 

written as 



\c n+ \ 2 = 4n 2 d 2 (- av lW+l&y 



= ^ 2 d 2 ( a u ii% + m, 



r(5) 



(45) 



with the coefficients l\f^ and 1^ given by Eq. i|tjt)[l . 
Using the dimensionless units, one has the radius R = 1, 
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the coupling constant b — 2/L so , and the wave number 
k = 2ttR/X, where A is the electron wavelength. Hence 
the semiclassical range of parameters corresponds to k » 
1. 

The residual polarization calculated from Eq. I|43|) is 
shown in Fig. 8. The P z curves for k = 20, 30, and 40 
are very close to each other and merge into the dashed 
curve. This curve coincides with the residual polariza- 
tion obtained from the semiclassical simulations in the 
previous section (Fig. 6). For k = 5, 1, and 0.1, the 
curves are plotted in the dotted, solid, and dash-dotted 
curves, respectively. All the curves, as expected, have the 
common asymptotic value 1 in the case of the weak spin- 
orbit coupling L so — ► oo. In the opposite limit, L so — ► 0, 
the behavior of P z is nonanalytic and not much reveal- 
ing. The strong oscillations in this limit increase with 
smaller wave numbers and signal about the appearance 
of large quantum beats in P z (t). This regime of L so is 
not interesting from the practical point of view because 
it implies unphysically large values of an for the typi- 
cal dot radius R = bOOnm. In the practically important 
regime of L so > 1 we note an apparent dependence of 
P z on k at k < 5. This is a quantum effect which is 
not observed in our semiclassical simulations. In semi- 
classics the particle velocity determines the speed with 
which P z (t) approaches to the residual value P z , but 
not this value itself. 
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FIG. 8: (Color online) The residual spin polarization P z vs 
L so with k = 20, 30, 40, (dashed), k = 5 (dotted), k = 1 
(solid), and k = 0.1 (dash-dotted). The dashed curve coin- 
cides with the curve from Fig. 6. 



V DISCUSSION 

Summarizing the above results of the semiclassical 
Monte Carlo simulations and quantum mechanical cal- 
culations we can draw the following picture of the spin 



evolution in semiclassical quantum dots. In the dots with 
regular classical dynamics the spin polarization does not 
decay to zero at long time and its residual value coincides 
with the quantum mechanical spin polarization averaged 
over an infinitely long time interval. At least, we were 
able to check such a coincidence for the circular dot. On 
the other hand, in dots with chaotic or random dynamics 
the spin polarization relaxes to zero with the relaxation 
time much larger than the DP relaxation time in un- 
bounded quantum wells. Such a decay down to zero can 
not be understood from the general quantum mechani- 
cal expression in Eq. (|32() . because it implies that the 
average of P z (t) over an asymptotically long time inter- 
val is zero. However, Eq. I|32l) predicts that this average 
is given by the first term in Eq. (|32[) . which is nonzero 
in general. Obviously, this contradiction is associated 
with quantum mechanical effects, which indicates that 
the semiclassical approximation is insufficient for anal- 
ysis of the long time polarization evolution. In disor- 
dered mesoscopic systems the statistics of their energy 
spectrum together with quantum interference effects give 
rise to the so called quantum dynamical echo [l5j which 
can contribute to the spin evolution at large times. This 
problem needs further study. 

The predicted spin evolution can be measured exper- 
imentally. For an InAs dot doped up to 10 n cm~ 2 , the 
time unit in Figs. 3-5 is about 1 ps if the dot size is 
L = 0.5/iTO. Hence, the spin polarization saturates to its 
residual value during first 20 ps and for L so = 1/im the 
difference in the long time evolution between chaotic and 
regular dots can be observed in the nanosecond range. In 
order to suppress all inelastic spin relaxation mechanisms 
EH3 El the measurement must be done at sufficiently 
low temperatures. The Rashba spin-orbit interaction can 
be strong in InAs based heterostructures, with L so down 
to several hundreds nm. Moreover, it can be tuned in a 
wide interval by varying the gate voltage 0]. 

In conclusion, we performed path integral semiclassi- 
cal simulations of spin evolution controlled by the Rashba 
spin-orbit interaction in quantum dots of various shapes. 
Our calculations revealed that the spin polarization dy- 
namics in QD's is quite different from the D'yakonov- 
Perel' spin relaxation in bulk 2D systems. Such a dis- 
tinction is not expected from the simple picture of the 
spin random walk, in particular when the rate of elec- 
tron elastic scattering on impurities in bulk is equal to 
the mean frequency of electron scattering from the dot 
boundaries. We have also found an important distinction 
between long time spin evolutions in classically chaotic 
and regular systems. In the former case the spin po- 
larization relaxes to zero within relaxation time much 
larger than the DP relaxation, while in the latter case it 
evolves to a time independent residual value. This value 
decreases with the stronger spin orbit interaction. We 
also analyzed the general quantum mechanical expression 
for the time dependent spin polarization. Using the exact 
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solutions of the Schrodinger equation with Rashba SOI 
for a circular dot, we calculated the average of the spin 
polarization over an infinitely long time interval and com- 
pared the result with the residual polarization from the 
Monte Carlo simulations. We found that the residual val- 
ues from these two approaches coincide, which confirms 
the results from the semiclassical simulations. On this 
basis, we conjecture that the nonzero residual value is a 
general property of regular systems. On the other hand, 
the spin relaxation down to zero in the Sinai billiard and 
circular dot with the rough boundary contradicts to what 
have to be expected from quantum mechanics. The long 
time memory due to the mesoscopic spin echo is assumed 
to be responsible for this contradiction. 



APPENDIX 

This appendix demonstrates a quantum mechanical 
calculation of the residual polarization P z , as it is defined 
in Eq. ll.'i.SII . The calculation of the exact eigenfunctions 
of the Hamiltonian in Eq. Q for the circular quantum 
dot can be found in Ref. |l6j . which is summarized in 
the following Eqs. I|46I50() . 

In order to calculate the residual polarization (|33|l , the 
wave function Vw(r) is expanded in the basis of the 
eigenfunctions given by Eqs. 
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We note that for 



a symmetric presentation, the functions f v and g v have 
different definitions from those in Ref. [16j. Inserting 
Eqs. (|35I36|) into the corresponding Schrodinger equa- 
tion we obtain the equation for f u and g v in terms of 
the dimensionless parameters £, e, and b defined in the 
previous section, 



[A„ + s] /„(£) - & V_(v +1 ) g u +i(£) = 0, 
[A v+1 +e]g v+1 (£)-bV +t ,f u (£)=0, 



with the Laplacian 



1 d 



and the nabla operator 

V±„ = ± 

The solutions (/</(£), <?i/(£)) 01 these equations are 
= di 



-o v Jv(*+0 + Mk-t) 



with the normalization constant d, the factors a v 
by Eq. i|38|) . and the wave vectors k± from Eq. 
These wave vectors obey the relations 

k^_k — £, k-\- k = 6, 



and k+ + fc_ = \/b 2 + 4e. 



(46) 



(47) 



(48) 



(49) 



given 
@2). 



(50) 



The quantized dimensionless energies e are determined 
by the zeros of the function in Eq. (|41fl . This function 
stems from the determinant of the equation system in Eq. 
(|4l)|l with the boundary conditions /„(£) = g v (Q = at 
£ = 1. Given a coupling constant 6, the n-th quantized 
value £„ with n = n(X, /i) is determined by the fi-th zero 
of Z v (e), where v and A are related by Eq. I|40|l . The 
allowed wave numbers k± are given by Eq. (|39|l with e = 
e n . They correspond to the two degenerate eigenstates 
of the n-th energy level. The first root of the function 
Z u (e) is zero for v = 1/2,3/2, ... and is a positive value 
for v = —1/2 (see Fig. 9). The larger the value of b, the 
larger the second root of Z v {e). 



0.6 
0.3 




v = -1/2 (solid), 1/2 (dashed), 3/2 (dotted) 




FIG. 9: The function Z u {e) for v = -1/2, 1/2, and 3/2. This 
function is singular at e = for v = —1/2. 

Substituting the wave functions in Eqs. (|35I36|I into 
Eq. Q33JI we obtain the residual polarization in the form 



p z = E / / [(K+l 2 -K-l 2 )(/, 2 (e)-^ + i(0) 

+2c* n „c n+ M0gv + i(0e l{2u+1)4 '} dc^Zdt. (51) 

For the initial wave function given by Eq. l|42f) . the 

constants c„+ can be expressed as 



Cn+ = 



<(r) \ 
^|*(r) e ^d 2 r 



ikr j2„ 



2tt 



i[k£t 



^-^-^Um^d^ (52) 



'o ^o 

where (f> and 9 stand for the angles of the vectors r and k 
with respect to the positive x-axis and k = |k|. After the 
shift of the angular variable from <f> — 9 to <f> the above 
integral transforms to 



C n + 



i{k£ cos(tj>)-u4 



e— - J /„(0£dW (53) 



Substituting t = 4> + 7r/2 and m = v into the integral 
representation of the Bessel function i 

Jm(z) = ^ ^ e ^ sia ^- mt Ut, (54) 
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2tt 



we obtain 

2ne l ^ /2 J„(z) = I eH'cnM-ufW^ ( 55 ) 
Jo 

By using this identity, Eq. I153H can be written as 

c n+ = 2ne 1 ^' 2 -^ [ Mk£)f (56) 
Jo 

By analogy, one has 

c n -=2ne^+^/ 2 -^ [ J„ +1 (k09„+imd£. (57) 



After integrating Eq. (|51|l over (j> and the direction 9 
of k (integration over 9 is similar to that over 6>k in Eq. 
(USD) the second and third terms in Eq. (|51|l vanish and 
only the first term remains. Introducing the parameters 

F*= I fimdZ and G n = [ <^ +1 (0^, (58) 
Jo Jo 

the final expression for the residual polarization can be 
written as 

pz = En(\Cn + \ 2 -\Cn-\ 2 ) (fn G n ) 

' E n (K + \ 2 + \c n -\ 2 )(F n + G n )- [ > 

For numerical calculations we explicitly wrote the norm 
of the normalized wave function ip(r, t) in the denomina- 
tor. In this form the expression in Eq. I|59(l is also valid 
for non-normalized wave functions, because the normal- 
ization constants d in the numerator and denominator 
are canceled with each other. 

The polarization P z in Eq. (|59ll is determined by the 
four integrals c n ±, F „. a nd G n . They can be calculated 
by using the formula [l£j 



UA^OM^O^ (60) 

_ I [nJ u (\l)J u +i(nl) - XJ u (kI)J v+ i(XI)] 

~ K 2 ~X 2 ' 

Consequently, the integrals in Eq. (|58|l can be written in 
the closed form 



F n =d z 
G„ = d 2 



all^-2a u I^+I^ 



a 2 I (1) +2a I {2) +I (3} 

a v 1 v+l ' La v 1 v+l ' 1 v+l 



(61) 
, (62) 



with 



IP 



MK) 2 - 


f J v +i{k+) 2 


uJ v (k + )J v+ i{k + ) 




2 


k + 


k-J v {k + 


J»+i{k-) - 


k + J v {k-)J v+ i{k + ) 




k 2 _ - 


k\ 


Ju{k-) 2 - 


f J,+i(k-) 2 


vJ u (k-)J v+ i(k-) 




2 


k_ 



By analogy, calculating the integrals in Eqs. (|56I57() 
we obtain 



4ir 2 d 



4tt z 



(_ ai/J (4) +7 (5)) ; (64) 

d 2 (a u li%+li%y, (65) 



with 



k + J v {k)J u+ x(k + ) - kJ v {k + )J u+ x(k) 



k\ 



k 2 



k^J„(k)J v+ i{k-) - kJ v {k^)J v+ i(k) 
k 2 -k 2 ' 



(66) 



For small b the spin polarization approaches to P z = 1 , 
as it must be in the absence of the spin-orbit interaction. 
It follows from the relation fc_ — k + = b w in Eq. I|5U|I , 
which results in /„(£) ~ 0, according to the definition 
in Eq. |@SJ|. Hence, the two quantities \c n+ \ and F n , 
which contain /„(£), vanish in P z . Therefore, the sums in 
the numerator and denominator of P z become the same, 
which gives rise to P z = 1. 

For large b, we have P z — > 0, which is due to the large 
difference between k + and fc_, namely, k- — fc+ = b 3> 1. 
According to the asymptotic behavior |l8| 



J v {x) = \ — cos 

TTX 



TT f 1 

X ~2 V + 2 



oi 1 - 

x 



(67) 



of the Bessel function at large x, the magnitude of the 
oscillating function J v (k + ) is much larger than J v (k_) 
by the order of \Jk- /fc+. Therefore, the leading terms 
of ii^ and 1^ in Eq. I|66[) behave like 



r(4) 



J v (k + )J u+1 (k) 
k 

Jy{k)J V +\{k-) 

k- 



(68) 



The first term is much larger than the second one. Con- 
sequently, both c n+ and c„_ are dominated by and 
have the same limit for large b. By analogy, F n and G n 
also have the same limit. Therefore, both |c n +| 2 — |c n -| 2 
and F n — G n in the numerator of Eq. I|59() become small 
and hence P z = for b — ► oo. 
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